Worse than a big rip? 
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We show that a generalised phantom Chaplygin gas can present a future singularity in a finite 
future cosmic time. Unlike the big rip singularity, this singularity happens for a finite scale factor, 
but like the big rip singularity, it would also take place at a finite future cosmic time. In addition, 
we define a dual of the generalised phantom Chaplygin gas which satisfies the null energy condition. 
Then, in a Randall-Sundrum 1 brane-world scenario, we show that the same kind of singularity at 
a finite scale factor arises for a brane filled with a dual of the generalised phantom Chaplygin gas. 
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INTRODUCTION 



In order to describe the current accelerated expansion 
of our Universe, use has usually been made of a fluid with 
equation of state p = wp, where p and p are the pressure 
and energy density, respectively, and w is & parameter 
which is expected to be fixed by observation. Observa- 
tional data [Jj suggest that w could be less than —1 and, 
in this case, the fluid is called phantom energy. It is well 
known that when this phantom energy is described as a 
K-essence field, from which the customary quintessential 
fluid is nothing but a particular case, the Universe will 
finish at the so-called big rip singularity Q, where the 
scale factor and the energy density go both to infinity 
and all the structures of the Universe are ripped apart, 
outside the horizon of each other. Phantom energy is 
most generally characterised by its violation of the clas- 
sical dominant energy condition. The big rip singularity 
can nevertheless be shown [3, |4( to be avoided when the 
universal fluid isgiven in terms of a generalised Chaply- 
gin gas (GCG) [5j, even when the dominant energy con- 
dition is violated, a case for which the equation of state 
takes on a more exotic, generally non-linear form. This 
conclusion can be extended to also encompass the case of 
the now dubbed dual GCG (DGCG) @. In this paper it 
will be seen that, even though the above conclusion can 
be maintained in the so-named dual version of phantom 
0], both GCG and DGCG models can predict another 
kind of future singularity, which might be called a big 
freeze singularity, as it takes place at a finite value of 
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the scale factor and infinite positive or negative energy 
density. Therefore, the big freeze singularity is different 
from the big rip singularity. In fact, the former singu- 
larity takes place at a finite scale factor while the latter 
occurs at an infinitely large radius. 

The content of this paper can be outlined as follows. 
In sec. II we shall consider how in the case of a GCG 
with positive energy density in a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) Universe the big freeze singu- 
larity (a singularity type III in the notation of Ref. [1] , see 
also Refs. [EllO|) can be obtained. The case for a DGCG 
occur ring in a brane world of the Randall-Sundrum I 
type [Tf|], with negative energy density is dealt with in 
sec. Ill, where it is shown that also in that case there 
will appear a future big freeze singularity. We finally con- 
clude and briefly comment on some particular aspects in 
sec. IV. 



II. PHANTOM GCG WITH A FUTURE 
SINGULARITY 

The phantomgeneralised Chaplygin gas was originally 
introduced in [13]. It satisfies the same equation of 
state as GCG M, i. e. 
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(i) 



where A is a positive constant and a is a parameter. If 
a = 1, Eq. |T]) corresponds to the equation of state of a 
Chaplygin gas. The conservation of the energy momen- 
tum tensor then implies 
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with B a constant parameter. It was noticed in Ref. [4| 
that if the parameter B is negative then the perfect fluid 
with the equation of state ([1]) cannot satisfy the null en- 
ergy condition; i.e. p + p < 0. Furthermore, it turns 
out that in this case the energy density grows, instead 
of red-shifting as the Universe expands. This kind of 
fluid was called a phantom generalised Chaplygin gas 
(PGCG). Moreover, the authors of 0] (see also [HI]) re- 
alised that if the PGCG parameter a is larger than —1 
then a FLRW Universe filled with this fluid would es- 
cape the big rip singularity because at large scale factors 
the Universe is asymptotically de Sitter. The authors of 
Ref. [H also pointed out that if a is smaller than —la 
FLRW Universe filled with this fluid has a maximum fi- 
nite size but did not work out the consequences in this 
case. 

Here, we will show that a PGCG (A positive and B 
negative in Eqs. ([1]) and ([2])) induces a singularity at a 



finite future cosmic time and for a finite scale factor in a 
homogeneous and isotropic Universe if a is smaller than 
— 1 . This singularity is different from a big rip singularity 
and we name it a big freeze singularity. In this case, it 
turns out that the PGCG energy density is an increasing 
function of the scale factor (see Eq. (J2J). In particular, p 
approaches at very small values of the scale factor 

(a — > 0) and blows up at a finite scale factor a max (see 
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Therefore, a FLRW Universe filled with a PGCG would 
face a future singularity at a finite radius [19|. Moreover, 
close to the singularity the cosmological evolution can be 
described as 
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For simplicity, we have assumed that the 3-dimensional 
spatial sections of the FLRW Universe are flat. Hereafter, 

2 

we shall use units such that ^ = 1, h = c = 1. In Eq. ^ 
tmax — t measures the cosmic time elapsed since the Uni- 
verse has a given scale factor a (a is close to a max ) till it 
has its maximum size; i.e. a = a max . It was previously 
pointed out in Refs [j, [l3| that a singularity at a finite 
scale factor can take place in a universe filled by a GCG. 
Here we have proven that this singularity not only takes 
place at a finite scale factor but also at a finite future 
cosmic time. 

On the other hand, a FLRW Universe filled with this 
fluid is asymptotically de Sitter in the past. In fact, 

a ~ a Q exp (A 2 V+*) tj , (5) 

where ao is a small scale factor. Moreover, the FLRW 
Universe starts its evolution at a past infinite cosmic time 
where a — > and p + p — + _ . In fact, the homogeneous 
and isotropic Universe super-accelerates; i.e. 

H = -\{p + p)>Q, (6) 

all the way until it hits the singularity at a — a max . In 
the previous equation, a dot stands for derivative respect 
to the cosmic time. We recall that the PGCG does not 
satisfy the null energy condition [J]; i.e. p + p < 0. 

In the next section, we will show that the same singu- 
larity; i.e. a big freeze, may show up in the brane-world 



scenario even for a fluid that satisfy the null energy con- 
dition. In fact, this can be the case for a dual of the 
phantom generalised Chaplygin gas. 

III. DUAL OF A PHANTOM GCG GAS AND 
RSI MODEL 

For a given phantom energy model described by a per- 
fect fluid, we define its dual as a perfect fluid which sat- 
isfies a similar equation of state with w < — 1 but fulfils 
the null energy condition 7] . Here w corresponds to the 
ratio between the pressure and the energy density of the 
perfect fluid. This definition of dual has not to be con- 
fused with the one used in Ref. [lOj . Despite the very 
weird property of the dual of a phantom energy model 
that its energy density is negative (2(| , it may supply an 
alternative to dark energy models in theories with mod- 
ified Friedmann equation like the brane-world scenario 
[l5| . In fact, this might be the case in Randall- Sundrum 
model with a unique brane (RSI) [111 ] where the modified 
Friedmann equation on the brane reads [l6| 

= />(' + £) • < 7 » 

For simplicity, we consider the flat chart of the brane. In 
the previous equation A is the brane tension which we 
assume positive. Consequently, the square of the Hubble 
parameter in RSI scenario is well defined for a negative 
energy density p as long as p < — 2A; i.e. as long as the ef- 
fective energy density of the brane is positive. Moreover, 
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the brane is super-inflating; i.e. the Hubble parameter is 
an increasing function of the brane cosmic time 



(8) 



It follows that a brane filled with the dual of a phantom 
energy is always accelerating; i.e. a > 0. The brane is 
super-inflating although the dual of the phantom energy 
satisfies the null energy density; ie. p + p > 0, because p 
is negative and smaller than — 2A. We remind the reader 
that the effective energy density of the brane is positive 
and is is well behaved (except at the singularity). This 
is a trademark of any dual of a phantom energy model. 

From now on, we will consider that the brane is filled 
with a dual of the phantom generalised Chaplygin gas 
(DGCG) [H]]. I n this setup, we will show how future 
singularities may show up at a finite scale factor and in 
a finite future cosmic time; i.e. a big freeze, on a brane 
filled by this fluid. Indeed, this is the case, for example, 
for a DGCG, therefore satisfying the polytropic equation 



of state ([T]) , whose energy density reads 



p= A 



B 
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where A is a negative constant, B a positive constant 
and n an integer. As can be seen, a commonness of the 
DGCG is that the parameter a is quantised; in this case 
l + a= -1/(1 + 2n). 

If n is positive, the brane faces a curvature singularity 
at a finite scale factor et max where 



(10) 



When the brane reaches its maximum size, the energy 
density diverges (cf. Eq. (J9j) ) . This feature is schemat- 
ically shown in Figs. [5] and [31 It can be easily checked 
that close to the singularity the scale factor can be ap- 
proximated by 
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Here t max is positive and constant. It corresponds to 
the future finite cosmic time (t max ) at which the brane 
hits the big freeze singularity. We would like to stress 
that the singularity is intrinsic to the brane. In fact, the 
bulk geometry is well-defined as it corresponds to two 
symmetric pieces of a 5-dimensional anti-de Sitter space- 
time split by the brane. 

On the other hand, depending on how large is the brane 
tension A (for positive n), the brane starts its evolution 
at a vanishing scale factor or at a finite non-vanishing 
scale factor ax- 
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In fact, if the initial energy density of the brane is smaller 
than — 2A; i.e. 2A < , the brane starts its evolution 

with a vanishing scale factor (sec Fig. In this case, 
the brane is asymptotically de Sitter in the past and its 



initial evolution can be approximated by 
a(t) ~ aiexjp(Hot), 
H - 



'L4|-(i 



|^4|-(l+2n) 
2A 



(13) 



where a\ is a very small scale factor and H Q is the ini- 
tial/primordial Hubble parameter. Therefore, the brane 
is born at an infinite past cosmic time because a van- 
ishes at t — > — oo. Now, if the initial energy density of the 
DGCG is larger than the brane tension; i.e. \A\ i+° < 2A, 
the brane starts its evolution with a non- vanishing radius 
a\ (see Fig. [3} . Therefore , the brane tension acts as a 
high energy cut off on p [22[ ■ This feature implies that A 
forbids an infinite past of the brane. Indeed, close to the 
initial scale factor a\, the cosmological evolution can be 
described as 



ax 
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where at t = t\, the brane radius coincides with a\. In summary, we have shown that a DGCG for n positive 
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(see Eq. [9]) induces a future singularity at a finite future 
cosmic time t max and for a finite radius of the brane a max . 
On the other hand, we have also shown that depending 
on how large is the brane tension, the brane starts its 
evolution with a vanishing radius at an infinite past or 
with a finite radius a\ at a finite past t\ (as measured by 
the cosmic time of the brane). At this respect, we would 
like to emphasise that we are considering only the phase 
where the brane expands. It is worth noticing that the 
expansion phase is preceded by a collapsing phase with 
negative Hubble parameter. It can be also seen that the 
latter phase starts with a past freeze singularity. 

Before concluding this section, we describe briefly the 
evolution of a brane filled by a DGCG whose energy den- 
sity is given by Eq. ((9|) for n negative (for more details 
see Ref. (6[). First of all, the minimum energy density 
of the DGCG must be smaller than — 2 A (see Fig. [4j; 
i.e. 2A < |A|~( 1+2n ) , otherwise the brane cannot have 
a Lorentzian evolution. Under this assumption, it turns 
out that the brane is asymptotically de Sitter in the fu- 
ture despite that the DGCG energy density is negative 
(see Fig. S]). At this respect, we recall that the modified 
Friedmann equation on the brane is quadratic on the to- 
tal energy density of the brane (cf. Eq. (JTJ) ) . In addition, 
the effective energy density of the brane is positive. The 
brane tension acts again as a high energy cut off on the 
DGCG energy density. Consequently, the brane starts its 
evolution at a non- vanishing scale factor a\ 1 previously 
introduced in Eq. (|12|) . Therefore, the brane tension ex- 
cludes those radii of the brane such that a m i n < a < a\, 
where the energy density vanishes at a m ; n . This cut off 
effect implies that the brane starts its evolution at a finite 
past cosmic time as measured by an observer confined on 
the brane. 



IV. DISCUSSION AND CONCLUSIONS 

We have seen how GCG and DGCG could generate a 
new big freeze singularity (a singularity type III in the 
notation of Ref. Q) which would be associated with a 
different way for the Universe to finish along a doomsday. 
At this doomsday the Universe would be infinitely full 
with phantom energy while it would have a finite size, 
in contrast with a doomsday a la big rip. Thus, nothing 
in the Universe could move around, producing a freezing 
for all the eternity. 

Close to a big rip singularity, the classical space-time 
breaks down and therefore a quantum analysis at large 
scale factors is required [13] • We expect similarly the 
quantum effects to be also important close to a big freeze 
singularity and consequently at intermediate values of 
the scale factor. In particular, the kind of big freeze 
singularity we have found in this paper needs an analysis 
of the GCG from a quantum point of view [l8( ■ We leave 
this interesting issue for a future work. 
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FIG. 1: Plot of p introduced in Eq. ([2} as a function of the 

scale factor. In this figure, we show the behaviour of p for 

1 + a < 0. A homogeneous and isotropic Universe filled with 

this fluid starts its evolution with constant energy density 
i 

A 1 + a at very small scale factor, therefore there is no big- 
bang singularity. The FLRW Universe super-accelerates all 
the way until it hits a singularity at a constant scale factor 
imax and in a finite future cosmic time. At the singularity 
the energy density diverges and consequently the Hubble pa- 
rameter blows up. 
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FIG. 2: Plot of p introduced in Eq. (|9} as a function of 
the scale factor. In this figure we show the behaviour of p 
(n in Eq. Q is positive) when the maximum energy den- 
sity is smaller than — 2A. Consequently, the brane starts its 
evolution from a vanishing scale factor. The brane super- 
accelerates all the way until it hits a curvature singularity at 
«max where the energy density blows up in a finite cosmic 

time tmax- 
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FIG. 3: Plot of p given in Eq. (|9} as a function of the scale 
factor. In this figure we show the behaviour of p (n in Eq. Q 
is positive) when the maximum energy density is larger than 
— 2A. Consequently, the brane starts its evolution from a non- 
vanishing scale factor a\ . The brane super-accelerates all the 
way until it hits a curvature singularity at a m&K where the 
energy density blows up in a finite cosmic time £ max - The 
dashed area corresponds to the region excluded from the evo- 
lution of the brane due to the cut off induced by the brane 
tension on the DGCG energy density. 




FIG. 4: Plot of p introduced in Eq. ((9} as a function of the 
scale factor. In this figure, we show the behaviour of p (n 
in Eq. ((9} is negative) when the minimum energy density is 
smaller than — 2A. Consequently, the brane starts its evolu- 
tion from a non- vanishing scale factor a\. The brane super- 
accelerates all the way until it becomes asymptotically de Sit- 
ter in an infinite future cosmic time (of the brane). 



